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PREFACE. 

This publication gives a general development of the 
theory of the Lambert conformal conic projection. It is 
intended to supplement the matter foimd in Special 
Publication No. 47 entitled, '*The Lambert Conformal 
Conic Projection with Two Standard Parallels." It is also 
supplementary in a way to Special Publication No. 49, 
which contains the Lambert projection tables for the region 
in France, and to Special Publication No. 62, which gives 
corresponding tables for the United States, since it gives 
as a whole the mathematical development of the theory 
upon which they depend. 

A short account of Lambert's life and work is given in 
the introductory paragraphs, followed by a few pages upon 
the subject of projections in general. 
2 
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GENERAL THEORY OF THE LAMBERT 
CONFORMAL CONIC PROJECTION. 



By Oscar S. Adams, 
Geodetic Computer, U. S, Coast and Geodetic Survey, 



Since this publication is to treat of some of Lambert's 
work, it is altogether fitting that it should be prefaced with 
a short account of his life, especially with a statement of his 
significance in the domain of map projections. Johann 
Heinrich Lambert was born at Mtilhausen in Alsace in 
1728. He was the son of a poor tailor and his education 
was entirely the product of his own exertions, due to a 
systematic course of reading. It was his regular custom to 
spend 17 hours pJer day in study and writing. At the early 
age of 16 he discovered, in computations for the comet of 
1744, the so-called Lambert's theorem. During the latter 
part of his life he resided in Berlin, where he was much 
honored for his abiHty. It was in the application of 
mathematical analysis to the practical problems of life 
that he especially excelled. His untimely death occurred 
in 1777 ia the forty-ninth year of his age. 

His contributions to mathematics were the series which 
bears his name, the conception of hyperbolic functions, his 
theorem in conies, and the demonstration of the incom- 
mensurability of TT. Both Lagrange and Gauss used parts 
of his work as points of departure for their investigations. 

Lambert's work in the field of map projections needs 
careful consideration. He was the first mathematician 
to make general investigations upon the subject of pro- 
jections. Those who preceded him in this work limited 
themselves to the development of a single method of pro- 
jection, principally the perspective, but Lambert con- 
sidered the problem of the representation of a sphere upon 
a plane from a higher standpoint and he stated certain 
general conditions that the representation was to fulfil, 

5 



6 U. S. COAST AND GEODETIC SURVEY. 

the most important of these being the preservation of 
angles or conformality and equal surface or equivalence. 
These two qualities, of course, can not be obtained in the 
same projection. 

Although Lambert did not fully develop the theory of 
these two methods of projection, yet he was the first to 
express clearly the ideas regarding them. The former, 
conformality, has become of the greatest importance to 
pure mathematics, but both of them are of exceeding im- 
portance to the cartographer. It is no more than just, 
therefore, to date the beginning of a new epoch in the 
science of map making from the appearance of Lambert's 
work. What he accomplished is of importance, not only 
for the generality of the ideas imderlying it, but also for 
his successful appUcation of them in methods of projection. 
The manner in which he attacks and solves any particular 
problem is very instructive. He has developed several 
methods of projection that are not only interesting but 
that are to-day in use among cartographers, the most im- 
portant of these being the conformal conic projection. 

The initial problem presented to us in map making is 
the representation upon a plane surface of the relative 
positions, sizes, and shapes of features that are foimd upon 
the curved surface of the earth. A perfect representation 
is impossible, since the surface of the earth being non- 
developable, can not be spread out in a plane. There are, 
however, many diflFerent ways of obtaining approximate 
representations, the theory and properties of which con- 
stitute the subject of map projections. 

The positions of the points upon the earth are usually 
defined by their latitude and longitude. Hence, if we can 
devise a suitable method of representing the meridians and 
parallels upon the sheet, the points can be plotted by their 
positions relative to these lines and the map can be con- 
structed. The term '* projection" is evidently used in a 
wider sense than that which is given to it in geometry. 
The majority of map projections are not projections in the 
geometrical sense — that is, perspective projections, orthog- 
onal projections, etc. — but merely a network of meridians 
and parallels that makes possible a one-to-one correspond- 
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ence between the places upon the earth and the points 
upon the map. 

Some of the things to be desired in a map are as follows : 

1. Preservation of the shapes of the countries. 

2. Preservation of the relative sizes of the countries in 
their representation upon the map. 

3. The distance between places should be in constant 
ratio to their distances as indicated upon the map. 

4. A great circle upon the earth should be represented 
by a straight line upon the map. 

5. The latitude and longitude of any place should be 
readily found from its position upon the map. 

6. The ease with which a projection can be constructed is 
also to be considered from the practical standpoint. 

Only part of these things can be attained by any given 
method of projection. 

The scale of a map in any given direction is the ratio 
which a short distance measured upon the map bears to 
the corresponding distance upon the siulface of the earth. 
The definition is limited to short distances because the scale 
of a map generally varies from point to point. It would 
of course be desirable that the scale of the map should be 
correct in every direction at every point and constant for 
all parts of the map. This is impossible, however, since if 
it were true, the map would be a perfect representation of 
the spheroidal siulface and could be fitted to it. In any 
given method of projection, therefore, some of the features 
to be desired must be sacrificed. 

The representation of the shapes of countries as nearly 
correct as possible is one of the most important functions 
of a map. A large country can not, of course, be repre- 
sented without some distortion when considered as a whole, 
but small areas may be mapped by similar figm-es. A pro- 
jection that preserves the similarity of small areas is called 
orthomorphic. The representation of one surface upon 
another so as to preserve similarity of elements has been 
called by mathematicians conformal representation. An 
orthomorphic projection is therefore a conformal repre- 
sentation of the spheroidal siulface of the earth upon a 
plane. Orthomorphic projections are in general not of 
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much use in map making unless the meridians are straight 
lines. 

The Lambert conformal conic projection fulfills this re- 
quirement, and it has lately been brought into especial 
notice by the fact that it is the projection that is used for 
the battle maps in France. In this projection, as is the 
case with any conic projection, the intersecting cone gives 
a better map than is given by the tangent cone. In the in- 
tersecting cone we have a shortening olf the scale between 
the standard parallels and a lengthening of the same at the 
top and bottom of the map. 

The following is a somewhat simplified development of 
the mathematical theory of the projection, including a 
determination of the necessary and sufiicient conditions 
for a conformal mapping. Use is made of the tangent cone 
only in order to determine the two limiting cases of the 
projection. 

It is well known that a plane curve can be expressed 
in parametric form in such a way that x = <t>(t) and y = 
^(0, </> and ^ being some known functions of the variable t. 
Thus the straight line passing through the point (a, b) with 
direction cosines a and P can be expressed in the form 

x = a + atyy = b-\-pt 

A circle can be given in the form 

• x^a cos tj y = a sin t. 

A space or skew curve can be similarly expressed as 

Such a curve may degenerate to a plane curve under cer- 
tain conditions; the obvious one is that one of the func- 
tions is identically zero. 
A straight line in space becomes 

x = a + atj y^b + Ptj z = c-{-yt 

A circular helix is given by the equations 

x = a cos tj y = a sint, z = bt 

When the coordinates are given in terms of functions 
of two variable parameters, the locus becomes a surface 
except under certain conditions that will be specified 
later. Thus, the tangent surface of a skew curve ex- 
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pressed by the parametric equations x = <t>{t), y = rp(t)j and 
z=x(f) is given by the equations 

the primes denoting differentiation with respect to t. By 
employing s and t as independent variables, the surface 
given is that swept out by a tangent to the curve as it 
moves along the curve. The surface of a sphere is given 
by the equations 

X = a sin ^ cos s, y = a sin ^ sin s, and z = a cos U 

To avoid some of the difficult questions connected with 
the theory of curves and surfaces, the functions that 
will be considered in this article will be limited to domains 
in which they are single valued, finite, continuous, and 
differentiable. These hmitations are permissible in the 
cases to be considered. Thus, in the equations of the 
sphere the variables are Umited to the domains 

— 2<^< +2 ^^^ —T^<s< +7r. 

In regard to the fimctions three cases may occur. In 
the first place all of the functions may be constant, thus 
giving a single point. In the second place, the functions 
may not be constant, but any two of them may be func- 
tions of the third; this would give a curve and not a surface 
since the three functions could then be expressed in terms 
of some single variable, ry, and the equations would become 

x = X(j])j y= r(r?), and z = Z(rj). 

As the third possibility, at least two of the functions are 
independent of each other. In this case the locus repre- 
sented is a surface. Of course any two may be inde- 
pendent, but all that is necessary is that at least two bear 
this relation to each other. 

Let and yp be the ones that are thus independent. 
This requires that their Jacobian is not identically zero, or 

du^ dv 



du dv 



4=0 
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It should be noted that a parametric representation of a 
surface is not unique. If such a representation has been 
found any number of other forms can be obtained by set- 
ting the u and v of the first form equal to two arbitrary 
independent functions of two other variables 8 and L 
Thus the equations may be u = X(s, t) and i;=m(«; t) with 
the condition 

+ 



d\ 


dX 


ds' 


bt 


bs' 


bt 



If one of the variables is held constant while the other 
varies, a curve is traced out upon the surface. Thus there 
is a set of curves for u = constant, and another set for 
v = constant. These curves are called the parametric 
curves and u and v are called the curvilinear coordinates. 
By a change of coordinates as indicated above a set of 
curves can be determined such that the curves, u = con- 
stant, are everywhere perpendicular to the curves, v = 
constant. In the equations of the sphere, s = constant 
gives a meridian great circle and t = constant gives a circle 
of latitude or a parallel. In this case the curves are per-, 
pendicular 'or orthogonal. In a more general way any 
functional relation between the variables u and v deter- 
mines a curve upon the surface. This can be expressed 
in the form of one coordinate as a function of the other, as 
v = g (u), or in the form F {u, i;)=0. As a further con- 
sideration any ordinary differential equation of the first 
order between u and v determines a set of oo^ curves upon 
the surface, since the integration of the same would give a 
functional relation between u and v containing one arbi- 
trary constant. 

The element of length upon a surface is given by the 

equation dS^=d7?+dy^+dz\ 

ou bv 

dy = ^du+^^dv 

ou ov 
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hence 



\ou ov ou ov ouovj 



+ 

To abbreviate the expression it is customary to use the 
symbols 

"du dv bu dv du bv 



H^y<^y<i^y 



With these the equation for the element of length becomes 

dS^^Edu^ + 2Fdudv-\-Gdv' 

If V is held constant, the element of length becomes 

dS^ = E du^, since dv = 

If the direction is taken positive in which u is increasing, 
we find 

dS^ = -yjE du 

Similarly for u = constant, dSy, = -yJO dv 

If V is constant dx = ^du, dy = ^dUj and dz = ^du. 

The direction cosines of the tangent to the curve v = constant 

dx 
are given by the equations a = -j^, etc. 

Thesebecome,a = -7= r— > P = —f=. r-^> and 7 =-7= r— 

' -y/E du -y/E bu -y/E du 

Similarly the direction cosines of the tangent to the curve 

u = constant are given in the form 

V6^ ov -^G bv -yJO bv 
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If 6 is the angle between these two tangents or the^ angl<> 
between the curves, we have from the known relation, 

cos = aa'+/3/3' + 77' 



7m^ 



dx dx dy by dz 




It follows, then, that the parametric curves u and v 
upon a surface on which F is identically zero form an or- 
thogonal system imless the -^^0 is identically zero at 
the same time. Since when F is zero neither E nor G 
is zero for the usual real surfaces, it follows that, when 

We shall now determine the conditions that are neces- 
sary and sufficient to permit a conformal representation 
of a surface upon a plane. One surface is said to be con- 
formally represented upon another when any infinitely 
small section of the first surface is represented by a similar 
figure upon the other. 

If a surface is to be represented upon a plane according 
to some law and x, y, z are the rectangular coordinates 
of the surface and u, v the rectangular coordinates of the 
plane, then for every pair of coor(finates u, v in the plane 
a set of coordinates x, y, z of the surface must correspond; 
that is X, y, and z must be functions of u and v. 

We may then consider that the surface with the para- 
meters u, V is to be represented on the plane with the 
rectangular coordinates u, v. The question then becomes 
under what conditions such a representation is conformal. 
Let us suppose that (u, v), (u+du, v-j-dv) and (u + 8Uf 
v-^-Sv) are three points upon the surface forming an in- 
finitesimal triangle. This will be represented upon the plane 
by the triangle (u, v), (u+du, v-^-dv) and (u + hb, v+Sv). 
Tliese two triangles must be similar but since their positions 
relative to each other do not come into consideration, 
no account is to be taken of signs either of sides or angles. 
Two things are then required, first that the sides of the 
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triangle extending from the point (u, v) of the surface 
must be in the same ratio as their representations in the* 
plane extending from the point (u, v) of the plane, and 
secondly the angle included between the sides of the tri- 
angle upon the surface must be equal to the angle included 
by their images in the plane. That is, the square of the 
elements of length upon the surface must be proportional 
to the square of the elements in the plane and the cosine 
of the angle between the lines upon the surface must 
equal the cosine of the corresponding angle in the plane. 

But ds' = Edv?-\-2 Fdudv^Odv", 

d8'-=E8u^ + 2 F bubv+Gbv". 

The quantities E, Fy and G would be the same in each case. 
They may be constants or functions of u and v but since 





Fig. 1. — Surface triangle and its representation upon a plane. 

each curve is infinitesimal and starts from the same point 
{Uy v), they would be constants for this point. On the 
other hand, in the plane, we have 

ds^ = du^'\-dv^ 

The first condition is therefore that the proportion 

Edu^+2 F du dv + G dv^ ^ du^ + dv^ 
E Su^-\-2 F 8u 8v-\-G dv^" M^-^dv^ 

may be identically true. This evidently requires that 
E=^G and F=0. Thus the parametric curves must be 
orthogonal. 

The cosine of the angle between the curves upon the 
surface must now be determined. The direction cosines 
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of the tangent at the point (u, v) of the curve dS upon the 
surface are determined by 

dx dy dz 
dS' dS' dS 

or dx J bx J 

>— au+5— dv 
ou ov 

etc. 



■y/Edu^ + 2 Fdudv^Gdv^' 

Also the direction cosines of the tangent at the point (u, 
v) to the curve hS upon the surface are in like manner 

bx^ hx ^ by ^ , by ^ 
ou ov ou ov 

^^/WW+TFWW+UW"' ^E 8u^ + 2 F 8u 8v + G 8v^' ^ ^' 
Consequently, if a is the angle between the curves, we have 

^^^^ ^(Edu^ + 2Fdudv + Gdv') (E8u^-]-2F8u8v + G8v'^ 

or 

E du au-f- F (du 8v-h8u dv)+G dv 8v 
cos « = 



^(Edu^ + 2 Fdudv + Gdv^) {E8u^ ^2 F 8u 8v + G 8v^) 

In the plane the corresponding direction cosines are 

/ du dv 

-yldu^ + dv^' -y/du^ -f- dv^ 

and 8u 8v 

'y/8u^ + 8v^' -yl^vF^^W 

If A is the angle between these two lines, we have 

. du8u-\-dv 8v 

cos A = —i= = 

^(du^ + dv^) (8u^ + 8v^) 

The second requirement is, therefore, that the equation 

Edu8u+F (du 8v + 8u dv) +G dv 8v 

V(£'du2 + 2 F du dv + G dv") (E 8u^ + 2 F 8u 8v + G 8v^) 

du 8u-{-dv 8v 

~ ^(du^ + dv^) {8u^ + 8v') 
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should be identically true. This requires as before that 
E=Gy and F=0. Thus of the two conditions only one is 
necessary and one is also sufficient since each leads to the 
same results. It is thus seen that the linear element upon 
the surface must be such that it can be expressed in the 
form dS^ = m^ {dvJ^ + dv^) . 

If we have a surface with the linear element d/S^ = 
Vdu^ -f Vdv^ in which F is a function of v alone and 
,r7 is a function of u alone, it can be expressed in 
the required form. Let us put it in the form d/S^ = 

-Tr+ -y )' By changing the parameters, we may let 

6 = I -j^ and rj = I ~T^' The linear element then becomes 

dS'^UVide^-^dri'). 

When the linear element is in this form, the surface is 
said to be expressed in terms of an isothermal orthogonal 
system of parameters and the net of u, v curves is said to 
form an isothermal orthogonal net. The surface is divided 
by them into an oo^ group of infinitesimal squares. 

After the surface has been expressed in this manner 
in terms of isothermal orthogonal coordinates, the general 
conformal representation of the surface upon a plane can 
be determined. This at the same time determines the 
general conformal mapping of the u, v plane upon another 
plane with point to point correspondence. This general 
representation is given by the equation x-\-i'y=f (u+iv), 

in which i denotes as usual V~ 1 and/ {u -f iv) is any ana- 
lytic function of the complex variable (u-fiv). The 
element of length dS^ =^w?{du^ +dv^) is represented in the 
plane by ds^^dy? +dy^. Let the curve C be represented 
by the curve C^ in the plane. If d represents the angle 
which the curve C makes with the curve v = constant, we 
will have dSv^mdu and 

. dSv du 

cos 6= 



dS -y/du^-^-dv"' 
dv 



sin ^=-7-— — ^— 
-Jdv?-\-dv^ 
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If 6^ is the angle that C^ makes with the axis of X, we have 

dx 



. . dy 



e<^=cos O+i sin 6=^ 



or 6**^= 



3 

du-\-i dv 
-^Jdu^ -\-dv^ 

du-\-i dv 



V du+idv 
du—idv 



du—i dv 
Similarly, 

^i-rosd 4-isin^= dx+idy _ I dx-^-idy 
e--cos(?,+ism(/, ^^2^dy' ^dx-idy 

^^ ^uk^ d^-^\Ay 
dx—i dy 
Therefore, 

^^^^ V du+i dv dx — i dy 
du — i dv ' dx-\-i dy 

or ^<(^-*i) ^*!f±i^ dx-idy 
dx-^-idy du — %dv 

But sincea5+ij/=/(t^4-iv), so also will x—iy^^fiu—w) 

By differentiation, we have 

daj + idy = (du + i<Zv)/' (u +tt?) 
or 

and 

dx-idy ^ . 

du—idv ^ 

the primes denoting differentiation with respect to the 
complex variable. By substituting these values, we have 

^ -f(u+iv)' 

If r and Fj are another pair of corresponding curves 
starting from the same point and their angles are denoted 
by <t> and ^j, we should find that e^*^*"^>) equaled the 
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same expression. Since this expression is constant for a 
given point, it follows that 

or 

so that 

6 — (t) = 6i — (t>i» 

That is, the angle between the curves upon the surface is^ 
preserved in the representation on the plane. If we had 
made use of the functional relation x+iy= F (u—iv) we 

should find the relation e2'<^+^> = jt/. . v 

Therefore in the same way as above 

^ + ^i = + «i 
or 

6 — <^ == <^i — wj , 

That is, the angles are equal but directed in the inverse 
sense, a diflference that is not taken into account in con- 
formal mapping. The most general conformal mapping 
is given then by the relations 

x-\-iy^f{u+iv) 
and 

x^iy:=F(u—iv) 

The two other forms are merely conjugates of these, i. e.,, 

x-'iy=f(u—iv) 
and 

x—iy^ F(u-\-iv). 

It is easily shown that these functional relations also^ 
preserve the proportionality of lengths. By differentiating' 
the relation a;+iy=/ (u+iv) we find that dx-\-idy=' 
(du+idv) f'(u+iv). By differentiating the conjugate of 
this, we obtain dx — idy = (du — idv) f (u — iv) . 
When these results are multiplied, we obtain 

(dxy - (idyy = L (duy - (idvy]f (u + iv)f (u - iv) 
or 

dx^ + dy^ = (du^ + dv^)f (u + iv)f (u - iv) 

But the element of length in the plane has the form 

ds^ = dx^ + dy^ 

67341"— 18 8 
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Also on the surface dS^ = m* (du^ + d^) 
By combining these relations, we find 

dS' 



or 



d^-^^f(u+iv)f(u--iv) 



ds ^ V/^ (u + iv)f' ju — iv) 
1[S~' m 



The two functions /' (u + iv) and /'(u— it?) are conjugate 
>complex functions and hence their product is real since 




Fio. 2.— Surface of revolution. 

both the sum and product of conjugate complex functions 
are real. The fact that the angle between the curves is 
held would in itself be sufficient to establish the comfor- 
mality of the representation, but we shall have need of 
the above expression which is called the magnification of 
the representation. 

If we have the equation of a plane curve given in the 
plane in the form x^u, y = 0, z=^f(u)j a surface of revo- 
lution with this curve as a meridian is given by the 
equations x=u cos v,y^u sin v, z =^f{u) . For this surface 
£=1 -l-[/(ti)p, F=0, and G^W, the prime denoting dif- 
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lerentiation with respect to u. The square of the element 
of length is therefore dS* = { 1 + [f'(u)y}du^+u^d'Vi' 



or 

u 



dS' = u2[ l + l/^'Wr ^^2 ^ ^^ 



We can let ^ ^ f Vl + LfWr 



=/ 



Then d/S^^u^Cdt^^'+dv^). The surface is thus expressed in 
terms of an isothermal orthogonal set of coordinates. 

Another test for a conformal mapping is that the func- 
tions must satisfy the Cauchy-Riemann partial differential 
equations. We have 

x + iy=f{u+iv) 

^^f(^u+iv) or ^-f(^+iv) 

Hence ^_^ j . ^__^, 

du~dv ^^ dv" du* 

With the form x + iy = F(u — iv) 
we obtain ^=;r'(„_i^,) 

Therefore ^= "I? *^^ ^=ll' 



20 TJ. 8. CX>AST AND GEODETIC SURVEY. 

These are the equations that determine that a given 
complex function is a function of a complex variable in 
the accepted sense. The first set preserves the angle 
between two given arcs in magnitude and sign. The second 
set results in reversing the angle, or the original angle and 
the image have opposite signs. Both x and y are functions 
that satisfy Laplace's equation for two dimensions. Thus 
with the first set 

b^x_ d^y 



bv? 


dvdtt 




dudv 


^y _ 


b'y 



But 

Therefore by addition 

Also d^y ^^y — n 

The same results may be found from the second set. 
From these results we may conclude that the most general 
conformal representation of one surface upon a plane is 
given by setting the complex variable in the plane equal 
to any analytic function of the complex variable formed 
from the isothermal orthogonal coordinates of the surface, 
or to any analytic function of the conjugate of this com- 
plex variable. In the first case direct equality of angles 
is found; in the second, the angles are equal but turned 
in opposite directions. 

If the surface to be represented is an ellipsoid of revolu- 
tion, the parametric equations may be chosen in the 
following form: 

a! = a cos M sin u, y = a sin M sin u, z = b cos u. 

a is the semi-major axis, h is the semi-minor axis, Mis 
the longitude and u the complement of the eccentric angle 
of the generating ellipse or the complement of the reduced 
latitude. 
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The element of length upon the spheroid is given by the 
equation dS^=dx^ -\-dy^ -\-d^. 
But 

dx = a cos JIf cos u du—asin Msm u dM; 

dy = a sin M cos u du -\-a cos M sin u dM; 

dz= —6 sin u du; 

hence the element of length becomes dS^ = a^ sin^ u dM^ -f 
(a* cos* u-\-b^ sin* u) du\ Thus E= a^ sin* u, F= o, and G = 

ct^ COS * u + 6* sin* te. If — ,— is put equal to c* the equation 
becomes 



(Z/S* = a? sin* -u [dif * + (cot* u -f 1 - c*) du*]. 

The equation of the generating elUpse is given by 

x = asin u; 
z = h cos u. 




Fia. 3.— Oenerating ellipse with angles L, y, and u. 

The tangent of the colatitude p, is the cotangent of the 

angle which the normal makes with the X axis. The 

tangent of the angle which the tangent makes with the X 

, , dz 
axis 13^- 

The tangent of the angle which the normal makes with 
the X axis, or the tangent of the latitude L is equal to 
dx 



dz 
But 



dx r ^ 

— -y- = tan L = x cot u 
dz 
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or 

tan L = cot V = t cot u 

hence 

- tan u = tan p 
a ^ 

or Vl — c^ tan u = tan p. 

, (1 — €2)cos2« . 1 

cos^ u = -^^ r — r^ > smce cos^ u =T~n — 5" 

1 — c^ cos^ 2>' 1 + tan'tA 

. , sin^p 

1 — c^ cos^ p 

Jul?, -in ^ 

^ cos^ u cos^ p 

or 

1 — e^ cos^ p 

(cot^ t^ + 1 - €2)du2 = jz — T--^^-^-r- 
^ ^ (1 — e^ cos^ py sm* p 

1 — c^ cos^ PL (1 "" ^^ cos^ p)2 sin^ j>J 

Let 

(1 — e^ cos^ 2>)sin 2> 
then 

r (?p € T c sin p (?p e f — 6 sm p <Zp 

"Jsinp 2J1 — ccosp 2 J l+€Cosp* 
By integration 

^ = logtan^ -Q ^^S (1— €COs p) + 2 ^^8 (l + €COS p)+log6^. 

If the Umits of integration are so chosen that G (the con- 
stant of integration) becomes equal to imity, we have 

r P /l + €COsp\ll 
^=logLtan2- Vi-ecosp^ J 
or 

c 



hence 



e' 



= tan2. / 1+^cosp y 
2 \1 — ccosp/ 
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The element of length now becomes 

1 — c^ cos' p\ / 

The parameters have been reduced to an isothermal 
orthogonal system. 

We can now determme any number of conf ormal repre- 
sentations of the spheroid upon the plane. All that is 
necessary is to make use of the relation 

x-\-iy^f (M±ie) 

f being any analytic function and either combination of 

signs being used. 

With the relation x^-iy^f (M-id), let/(v) = £^«». 

We now have 

x+^t/=Ze ^L^wij vi-fcosp/ J 



^ii 



= ^>-[-f (1^^)^] 



XBiiM 



= Ztan' 2. / I + ^cosA Vcos IM+i sin Im). 

2 \1— €COSp/ \ / 

By equating the real parts and the imaginary parts, we 
obtain 

«=Ztaii'|. ( l+'^^P )^ cos IM 
2 \1 — c cos p/ 

y=Ztan'|. (]±lSS^)Km I Ml 
^ 2 \1 — c cos p/ 

In this projection the parallels become concentric cir- 
cles. The equation for the radius is 



^ ^ 2 \1 — € COS p/ 



h 

2 



The meridians are represented by radii of these concentric 
circles. This method of projection is the one known as 
Lambert's conformal conic projection, first developed by 
John Heinrich Lambert, in his ''Beitrage zum Gebrauche 
der Mathematik," Berlin, 1772. It was later fully dis- 
cussed by Gauss. 
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If an angle z is assumed such that 

tan|=tan£. ('i+l^^V 
2 2 VI — « cos f) 



this angle will be very nearly equal to the complement of 
the geocentric latitude. 




FlO. 4.-»0eneratiiig ellipse with angles L and L', 

From the equation of the generating ellipse (p. 21), the 
tangent of the geocentric latitude is given by 



But 



or 



hence 



- =- cot u = tan L' 

X a 



tan i=T cottt 



cot u==— tan Z, 
a 

tan i' = - cot tfc = -a tan X, 



a and h being the semimajor axis and semiminor axis, 
respectively. 

Then to a sufficient degree of approximation z — ^—L\ 

The value of z can be computed rigidly very conven- 
iently by assuming an angle q such that cos q==€ cos p. 
Then, since 

2-"Vl"- 



cot 



cos q 
cos q 



tan^ — tan ^ cot* ^ 
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or 

log tan K = log tan ^+€ log cot |* 

However, the approximate formula determines z to within 
a few tenths of a second. By using this auxiliary angle 
the equations become- 

x = -K" tan^ 2 ^^s ^-^ 
t/=Z'tan'2 sin IM 

r=Z'tan^2- 

With these values x is reckoned downward from the center 
of the concentric circles and y to the right of the central 
meridian if Jf is reckoned positive in that direction. 

K and I are as yet left arbitrary constants. I may be so 
determined that the ratio of the lengths of any two arcs of 
parallels on the map may"" be equal to the ratio of the 
lengths of the arcs that they represent. If N is the radius 
of curvature perpendicular to the meridian, or the length 
of the normal prolonged to the minor axis, a radian of the 
parallel L^ has the length iVj cos L^\ in the same way the 
length of a radian of parallel L^ is iVj cos L^* Conse- 
quently, the ratio of the lengths of these two arcs is repre- 
sented by 

N^ cos L^ 

N2 cos L2' 

Since the A factor^ in the tables for geodetic positions is 

equal to 

1 



the ratio becomes 



iVsinl" 
A2 cos Li 



A^ cos L2 

The arc upon the map that represents the radian of parallel 

z 
ii has the length Ir^^lK tan^ ^. The radian of parallel i, 

1 See United States Coast and Geodetic Survey Special Publication No. 8, entitled 
•'Formulae and Tables for the Computation of Geodetic Positions.'' 
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is likewise represented by Ir^^lK tan^ -^' The ratio of 
lengths will be preserved if we have 

tan ^\^ 

tan I /"Acoai, 

, _ log cos Zi — log cos L2 — log ^t +log A2 

log tan Y - log tan —^ 

K may now be determined so as to hold not merely the 
ratio of the arcs of parallels L^ and ij ^^^ ^^so to hold the 
exact length of these parallels. This is an excellent method 
of determiilation for mapping such an area as that of the 
United States. In this way we should have 

IK tan^ ?r = iVi cos ii= -1 — = — hyj 
2 ^ ^ ^ism V 

hence r t 

j^ cos L^ cos L2 



^2 



A, sin V I tan' ^ A^ sin 1" I tan' ^ 

The twofold determination serves as a check on the com- 
putation. 

With this determination of I and K, we shall compute the 
expression for the magnification at any point. We 
employed as the form of the function / 

From the equation of the linear element on the ellipsoid 

on page 23 

^ ^ a sm p 

m= , — ^ 
yl — e^ cos^ p 

a sm p 
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But 

2 \1— € COS p/ 

therefore 

Kl tan^ I 



fc = 



(T sm p 



1/1 +€ cos 2)\^ /:; r r— 

( z. ]^Jl—e^ CO^^ p 

\1— € COS 2>/ 



or 



p^ iV^ COS Ln COS i^ 

r^ being the radius of the circle on the map that represents 
the parallel of L^ and p„ being the radius of the parallel. 
The last form is obvious from the conditions. 

If the parallels to be held are chosen about 3^ of the 
distance from the bottom and the top of the area to be 
mapped, the proper balance will be preserved. The 
upper and lower part of the map will then be about as 
much too large in scale as the central part is too small. 
The scale along L^ and L^ will be exactly correct. With 
this value of K one can tell how much any parallel is in 
error of scale by computing a radian of the parallel and 
the length of the arc which represents it on the map. This 
is just a statement of the equation 

Pn 

With this projection a map could be made of an area 
such as that of the United States so that it would not be 
in error of scale in any part of it by more than 1 1/5 per 
cent.. A poly conic projection of the same area is in error 
of scale by as much as 6^^ per cent in some parts. A 
Lambert projection for the United States to be evenly 
balanced should hold parallels 29*^ and 45*^. The scale 
would then be just about 1 per cent short along the 37*^ 
parallel and 1 1/5 per cent long along the 49° parallel. 

In Special Publication No. 52 of the United States 
Coast and Geodetic Survey, are given tables of coordinates 
for a Lambert projection of the United States. The 
standard parallels were chosen as 33*^ and 45° to lessen 
the scale distortion of the middle section of the country. 
This scheme gives about 23^ per cent scale distortion 
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along the parallel of 25° as the greatest on the whole map. 
Only a small part of Florida and of Texas are affected by 
this increase of scale error. 

The coordinates for mapping the parallels are most con- 
veniently computed using as origin the point where the 
parallel crosses the central meridian, the central meridian 

being the F axis and a perpendicular to it 
the X axis. All the formulas required for 
computation are as follows: 

tan i' = — i tan L 
a? 

log tan -^ —log tan -^ 

cos L^ cos Z>2 

A, sin V I tan' -| A^smVl tan' ^ 




Z= 



x = r sin IM 
no. 7.— Representa. J/ = r (1 — COS LM) = 2r sm^ -9- = a; tan -«- • 

tion of a small geo. ^ ^ 

detic trapezoid. t/ -x • j • j x 1 xt_ x x- 

if it IS desu-ed to make the computation 
considering the earth as a sphere, it is only necessary to 

IT 

let € = 0. z then equals the polar distance or z=p=^ — L; 

1 



A^=A^ 



When these values are inserted in the 



a sin V 
above formulas, the correct forms are given for the sphere. 

The difference of the radii gives the spacing upon the 
central meridian. If the parallel at the top and the one 
at the bottom are constructed by determining the coordi- 
nates of their intersections with the meridians, the me- 
ridians can then be drawn. They can then be subdivided 
as was done in the case of the central meridian. In this 
way the coordinates of the other parallels can be deter- 
mined without computation. 
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If the two parallels to be held approach each other indefi- 
nitely, we shall have to determine the hmiting value of I. 
This can be done in a number of diflFerent ways, and 
in fact the computation was made in six diflFerent ways 
with the same result in each case. The value that was 
obtained is the same as that obtained by Gauss and cited 
by Forsyth in his "Theory of Functions of a Complex 
Variable.'' Especial pains were taken with this, because 
Germain in his "Trait6 des Projections'' obtains an 
erroneous result. The value is here determined by the 
same method that Germain said that he used, not because 
it is the simplest but to illustrate the fact that the correct 
result can be so obtained. 

The problem is to so determine I that CD may have the 
same ratio to BD that AB has, or to determine the place 
where two successive parallels are held at the same ratio. 
This is the same as saying that the magnification shall be 
a minimum at this point and can be true only at this point. 

length of BD on the map is equal to dr 

length of BD on the earth = LLE — _ 

(1— €^ cos' p)"2 

length of CD on the map = Z(r+dr)dif 

length of CD on the earth 



-[ 



L V 1 — € cos' p \ yi — e' cos' p/S 

asm p a cos p dp a^ sin' p cos pdp' y 

+ . - , 3 ^M 

yl — €'cos'2> ^\ — ^ 0,0^ p (1 — €' cos' j>) 2 J 



=[ 



a si n p a(l — «') cos pdp^ 

V 1 — €' cos' p (1 — €' cos' p)2 J 



To meet the required conditions, we must have the propor- 
tion 



t g siny a(l — €') cos ydp ") 



• Z(r+dr)dif ^L yi -6' cos' p (1 -c' cos' j>)Tj 
dr a(l-€^)dp 



(1 — €' cos' p) 
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Therefore, since 

— (1 — t*) dp 

Ir'^il— t' cos' p) sin p' 

l = coa p. 
The magnification is least at 

p ™ COS"' Z or at i = sin""' I. 
If this least value of the magnification is denoted by fr, 
we have 

If Z is taken equal to cos p, and ^is detennined so as to 
hold the length of the parallel of colatitude p, we shall 
have the case of a tangent cone. 



FiQ. S.— St«reagTaphic pcoJectloD of ths northern hemisplien. 

Germain obtained the erroneous value 

J cos p — ^ COS p(l — 3 sin' p) 
*= 1-*' 

If I becomes equal to unity, we have po = o and the 
tangent cone becomes a tangent plane at the north pole* 
The equations now become 
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r=Ztan|- 

x=r cos M 
y^r sin M, 

This gives a projection for the spheroid analogous to the 
stereographic projection for the sphere. With the sphere 
the value of € is zero and z becomes the polar, distance. 
We then have a perspective projection from the south 
pole upon the plane tangent at the north pole. 

If I becomes equal to zero the tangent cone becomes a 

cylinder tangent at the equator, and po equals -^. To 

Jul 

determine the values it becomes necessary to evaluate 
some limits. K becomes infinite but so that Kl is finite 
and K—x finite: 

iVo cos io 



Z= 



ZZ = 



When 



Po — 2o — 2 



/ tan^ § 

Ad 

iVp cos Lq 
tan' I 



TT 



then 

lim Kl=No=a 

In the general formulas 

2 \1 — € COS p/ 

y vu^ 2 \1 — 6 COS pj 

To evaluate the limit, let us write x in this form 

or on developing the exponential 

*-iroo8ZM{l + Zlog[taii|. (l±l^)i]+0(P) 
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0(jP) denoting terms with P as the lowest power of I. But 

hence 

Xo-x=K(l-cm IM) - Kl cos Of log [tan |- ('l±i^£ ^Vl 

-KcoslMOQ?). 

TT 

Taking the limit when 2>o == k » l=='0 



X, 



-x = alog cot-^. / ^ <^Qsp \ 

L 2 \H-€ cos p/ J 



Denoting Xq — x by x and substituting for p its value ^ — Lj 

we have 

x = alog„[tan(^+j). ([^l||^)i] 

logrt denoting the Naperian logarithm. 

y=m tan^ f • ('I +i£2LE)^ ^JH^. 
^ 2 \l-€COsp/ Z 

Taking the limit when 7>o==k> Z=o 

M being of course expressed in arc. Interchanging x and 
y to give the coordinates as usually plotted, we have 

x^aM 

y = alog„[tan(^+j). ([^llg^)^]. 

This is a projection of the spheroid analogous to the 
Mercator projection for the sphere. If € becomes zero we 
have the Mercator projection of the sphere. 

We thus find that the stereographic and the Mercator 
are special cases of the Lambert projection and are there- 
fore conformal. 

Certain points in a conformal projection may be singular 
points at which the conformality fails. This is the con- 
dition at the pole in the Lambert conformal conic projec- 
tion. The angles between the meridians are not pre- 
served, since the angle between two meridians is ZJf instead 
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of M as Upon the earth. At such points, if w—f(z) is the 
complex relation, -^ is equal to zero or to infinity. The con- 
formality fails at the poles in a Mercator projection but is 
preserved at the central pole in a stereographic polar pro- 
jection. Since the angle at the center of the system of 
circles in the Lambert projection is equal to IM, the 360° 
in longitude is mapped on a sector of a circle with central 
angle equal to 360° X I. • Since Z is in the usual cases less 
than unity, the central Angle will be less than 360°. If I 



Fio. B.— MeraalOT prcijeolion. 

is equal to % , the angle of the sector would be 270° ; for I 
equal to ^ the angle would be 240°. 

In the war zone in France, the maps are made in the fol- 
lowing mamier. Instead of using the exact formulas that 
have been developed, approximations are employed. A 
cone tangent at 65 grades (49° 30') is first determined with 
the parallel of 55* as the central parallel of the map. 
Along the central meridian the parallels are spaced from 
the formula 

6po' 
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P being the distance along the meridian on the earth and po 
being the mean radius of curvature at 55^. This formula 
for Ar is a Taylor series development of Ar from the rigid 
formula for r, correct to the third power of p} 

The radius for the parallel of 55* is taken a^ No cot 55*, 
No being the radius of curvature perpendicular to the 
meridian. The values of Ar being added to or subtracted 
from this radius give the radii of the other parallels. A Jf 
sin 55* gives the arc along the parallel correspondiag to 




Fig. 10.— 'Graticule for the Lambert projection of the northern hemisphere with i— J^. 

the difference of longitude AM reckoned from the central 
meridian. • 

After these values are computed the whole scale is re- 
duced by 1 part in 2037. This gives us approximately a 
cone intersecting at the parallels of 53* and 57*. The 
whole map is then covered with a system of kilometer 
squares with origin at latitude 55* and longitude 6* east 

1 See United States Coast and Geodetic Survey Special Publication No. 47, p. 13. 
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of Paris. The lines north and south are all parallel to the 
meridian of 6*^ and the east and west lines are perpendicular 
to the same. A great circle for a limited region is a straight 
line within the lunits of scalmg. Since the map is con- 




FiQ. 11. — Graticule for the Lambert projection of the northern hemispheae with Z>-^. 

formal, a chart made upon this projection is of great use 
in determining the direction for gunfire. The scale is also 
preserved constant within the error of scaling, so that the 
projection is excellent for the determination of both direc- 
tion and distances. 
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